Novel universality classes of coupled driven diffusive systems 



CO 

o 
o 

CN 
O 

Q 



Abhik Basu 1,2 and Erwin Frey 2 ' 3 
1 Poornaprajna Institute of Scientific Research, Bangalore, India, 
2 Abteilung Theorie, Hahn-Meitner-Institut, Glienicker Strasse 100, D-14109 Berlin, Germany, 
3 Fachbereich Physik, Freie Universitat Berlin, Arnimalle 14, D-14195 Berlin, Germany 

(Dated: February 2, 2008) 

Motivated by the phenomenologies of dynamic roughening of strings in random media and mag- 
netohydrodynamics, we examine the universal properties of driven diffusive system with coupled 
fields. We demonstrate that cross-correlations between the fields lead to amplitude-ratios and scal- 
ing exponents varying continuosly with the strength of these cross-correlations. The implications of 
these results for experimentally relevant systems are discussed. 

PACS numbers: 0.5.70.Ln, 64.60.Ak,0.5.40.-a 



c3 



-a 

G 

o 
o 



> 



CN 
CO 

o 

-I— > 



i 

C 

O 

o 



X 



Recently significant advances have been made in clas- 
sifying the physics of non-equilibrium systems at long 
time and length scales into universality classes. It has 
been shown that standard universality classes in criti- 
cal dynamics are quite robust to detailed-balance vio- 
lating perturbations Q. Novel features are found only 
for models with conserved order parameter and spatially 
anisotropic noise correlations. In contrast, truly non- 
equilibrium dynamic phenomena, whose steady state can 
not be described in terms of a Gibbsian distribution, are 
found to be rather sensitive to all kinds of perturbations. 
Prominent examples are driven diffusive systems |2| and 
diffusion- limited reactions For example, one finds 
that for the Kardar-Parisi-Zhang equation anisotropic 
perturbations are relevant in d > 2 spatial dimensions, 
leading to rich phenomena that include novel universality 
classes and the possibility of first-order phase transitions 
and multicritical behavior p|. 

In this letter we study driven non-equilibrium pro- 
cesses described by a set of dynamic variables whose 
dynamics is given in terms of coupled Langevin equa- 
tions. Prominent examples include the dynamic rough- 
ening of strings movingin random media |5| , sedimenting 
colloidal suspensions p and crystals 0, and magneto- 
hydrodynamics (MHD) 8]. Our goal is to investigate 
and elucidate some of the dramatic effects of symme- 
tries of correlation functions on the universal proper- 
ties of such systems. We focus on models with two 
vector fields, u(x, t) and b(x, t), as hydrodyamic vari- 
ables. The quantities of interest are the two auto- 
correlation functions, C,"(x, t) = (itj(x, t)uj(0, 0)) and 
C|-(x,t) = (fc(x,t)& 3 -(0,0)), and the cross-correlation 
function C£j(x,t) = (uj(x, t)bj(0, 0)}; indices i, j refer to 
cartesian coordinates. All these quantities are tensors, 
whose symmetry properties depend on the model under 
consideration. We are interested in systems with transla- 
tional and rotational symmetry, and inversion symmetry 
such that u is a polar and b is an axial vector. 

In the first part of the letter, we will consider 
a one-dimensional Burgers-like model || of magneto- 
hydrodynamics and its c£-dimensional generalization [l(Tj | 



dh 



^V 2 u + f, 



^ + A 3 V(u • b) = /iV 2 b 
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Here \ are coupling constants, v and are the dissipa- 
tion coefficients, and f and g are external stochastic forc- 
ing functions. These equations are simplified versions of 
the dynamical equations governing the time evolution of 
the velocity u and the magnetic field b in a magnetized 
fluid (MHD). They are constructed in the same spirit 
as Burgers equation from the Navier-Stokes equation. In 
the second part of the letter we will discuss the advection 
of a passive vector b, where Ai = A2 = 0. The simplic- 
ity of such a model will allow us to explore higher order 
correlation functions. 

For Langevin equations describing processes relaxing 
towards a thermal equilibrium state the correlation func- 
tions for the noise have to obey detailed balance condi- 
tions. In non-equilibrium models there are no such re- 
strictions. As a minimal requirement one might ask that 
the noise terms f and g in the Langevin equations obey 
the same symmetries as the correlation functions for the 
hydrodynamic fields. Since u is a polar vector and b 
is an axial vector, (uj(k, t)v,j(— k, 0)), (&j(k, t)bj(— k, 0)) 
are real and even in k, but the cross-correlation func- 
tion C*j(k, t) — (ui(k, t)bj(— k, 0)) is imaginary and odd 
in k [8j. Then, assuming Gaussian distributed conserved 
noise with zero mean, the noise correlation functions have 
to be of the following form 

(/ l (k,<)/ J (-k,0)) = 2k i k j D^(k)5(t) (3) 
<Si(k,t) Sj -(-k,0)) = 2k l k J D^(k)S(t) (4) 
(Mk,t) 9j (-k,0)) = 2 t Df\k)S(t) (5) 

where the noise variances (k) are even and D * ^° (k) 
is odd in k, respectively. Equations (3) and (4) are invari- 
ant under inversion, rotation and exchange of i with j. 
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We take the noise cross-correlation, Eq. (5) to be invari- 
ant under inversion, but we allow it to break rotational 
invariance or symmetry with respect to an interchange of 
the cartesian indices i and j . 

We are interested in the physics at long time and length 
scales. Then all the correlation functions C(x,t) are ex- 
pected to obey scaling relations of the form 



C(x,i) =x 2x C(t/x z ). 



(6) 



Since we have two independent fields u and b there could 
in principle be two different roughness exponents Xu,b- 
Due to Galilean invariance, however, none of the non- 
linearities in the equations of motion rcnormalize, and 
one gets Xu = Xb = X = 2 - z H®. 

Symmetric Cross- Correlations. If both the fields 
are irrotational, one can introduce two scalar fields h and 
4> such that u = V/i and b = V<ji>; note that (j> is actually 
a pseudo-scalar. Then Eqs.(|ll2|l become identical to a 
model of Erta§ and Kardar JJj describing the dynamic 
roughening of directed lines, 



(7) 



the scaling form Eq. In Fourier space the scaling 
form reads for the self energy, S(As, u>) = Tk z a(w/k z ), 
and for the generalized kinetic coefficients Dh(k,u>) = 
D h k- d - 2 *d h {u/k% D^k,u) = D^k- d - 2 ^{u/k z ), 
£> x (k,w) = sgn(k)D x k- d - 2 ^d x {u/k z ). To solve this 
set of coupled integral equations we employ a small x~ 
expansion This requires matching of the self en- 

ergies and correlation functions at uj = 0. With the 
zero- frequency expressions E(k, 0) = Tk z , Dh(k, 0) = 
D h k- 2x - d , £>0(k,O) = D^ 2 ^, one finds for the one- 
loop self-energy (we take Ai = A2 = A3 = A without any 
loss of generality) 



D h A 2 



s d 1 

(2ir) d 2d 



D4 



(9) 



and for the one-loop correlation functions, 
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D^X 2 2 (2ir) d d-2 + 3 X 
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(10) 



^ + A 3 (V/i)(V^) 



(8) 



where f = V^/j and g = V??0. The cross-correlation 

function Df^ is now symmetric in the tensor indices 
and (/i(k, u)(f>(— k, 0)) is imaginary and odd in k. If, 
in addition, we require rotational invariance, the cross- 
correlation function would vanish. This is the case con- 
sidered in Ref.||. For a truley non-equilibrium model 
there is, however, no physical principle which would ex- 
clude a finite cross-correlation term a priori. Hence we 
allow for a non-zero (??/i(k, 0)77^ ( — k, 0)), which then ex- 
plicitly breaks rotational invariance, and explore its con- 
sequences for the dynamics. 

We have determined the roughness exponent x an d 
the dynamic exponent z employing a lowest order self- 
consistent mode coupling scheme and a one-loop dynamic 
renormalisation group calculation. Perturbation theory 
is formulated in terms of the response and correlation 
functions for the fields h and <fi. They are conveniently 
written in terms of self-energies u>) and generalized 
kinetic coefficients D(k,u>). For simplicity we assume 
that v = fi; in MHD this would correspond to a system 
with magnetic Prandtl number P m = fj,/v = 1. Then 
there is only one response function and it can be written 
as G^^(k, cj) = iu — S(k, u). Then, correlation func- 
tions are of the form, C a (k,uj) = 2D a (k,uj)\G(k,u})\ 2 
for a = h,(j) and C x (k,w) = 2iD x (k, w)\G(k, uj) \ 2 
for the cross-correlation function. In diagrammatic 
language lowest order mode-coupling theory is equiv- 
alent to a self-consistent one loop theory. The ensu- 
ing coupled set of integral equations is compatible with 



Here is the surface of a d-dimensional unit sphere. 
From Eas. i|l(Jfl we find 



£±) +2^^ + 11-1 = 0, 



D 



(11) 



where N = (D x /Dh) 2 defines an amplitude ratio. In 
the Eq. i|ll|) , the domain of N is determined by the range 
of real values for D^/Dh starting from 1 (for N = 0). 
Thus for small N we can expand around and look for 
solutions of the form D^/Dh = 1 + aN, such that for 
N = we recover Dh = D^ (the result of Ref. |5|]). We 
obtain a — —2, i.e., 



Dt/D h 



1-2N . 



(12) 



implying that within this approximate calculation TV can- 
not exceed 1/2, i.e., D x < Dh/2. An important con- 
sequence of this calculation is that the amplitude ratio 
Dfj, j Dh is no longer fixed to 1 but can vary continuously 
with the strength of the noise cross-correlation amplitude 
D x . These results are confirmed by a one- loop RG for 
the strong coupling fixed point in d = 1 . In addition, Eq. 
(12) is valid at the roughening transitions to lowest order 
in a d = 2 + e expansion. 

In contrast, the scaling exponents x an( i z are not af- 
fected by the presence cross-correlations. We get X = \ 
and z = I in d = 1 dimensions, x — ~0{e) 2 and 
z = 2 + 0{e) 2 at the roughening transitions in a d = 2 + e 



expansion, and x — 2 — I and z — 



I for the strong 



coupling phase. Note that the values for the strong cou- 
pling exponents are obtained within Bhattacharjee's 01 
small-x expansion, as described above. There is still an 
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ongoing debate whether those values for the exponents 
actually correspond to the KPZ strong coupling case (see 
Ref.|13| for a discussion). 

We have verified our analytical results in d = 1 by 
numerical simulations of both a coupled lattice model 
with cross-correlations, and direct numerical simulations 
of the model equations (|TJ and (0) . Our numerical results 
explicitly demonstrate the dependence of the amplitude- 
ratio on the cross-correlation function amplitude. Details 
will be published elsewhere [l3j |. 

Anti- Symmetric Cross- Correlations.- In the preced- 
ing paragraph we have restricted ourselves to irrota- 
tional fields. If the vector fields a = u, b have the 
form a = V x V a + VS a , with vectors V a being cross- 
correlated but scalars S a uncorrelated then the vari- 
ance D* satisfies D*(k) = -£>?J(-k) = £>£(-k) - 
— [D*(k)]*.. This is the antisymmetric part of the 
crosscorrelations. The noise strength D x is defined 
by D£(k)D£(-k) = D 2 X fc 4 . In the scaling limit, the 
self energy reads ui) = Tk z a(uj/k z ), the correla- 
tion functions are Cfj(k,uj) = kikjD u k~ d ~ 2x d u (ui/k z ), 
C^(k,uj) = kikjD b k- d - 2x d b {uj/k z ), and the anti- 
symmetric part of the crosscorrelation function reads 
C§(fc,w) = D?.(k)fc- 2 ^- d . 

Following methods used for the symmetric crosscorre- 
lations, we obtain the analogues of Eqs.© and 1)1 
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(14) 



1 at the fixed 



(D x /D h ) 2 . Hence no 



Equations and 1(11)1 give D u /D b 
point for arbitrary values of N = (D x /Uh 
restrictions on N arises from that. In contrast to the 
effects of the symmetric crosscorrelations, the exponents 
now depend continuously on N. To leading order, we get 
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These exponents presumably describe the rough phase 
above d > 2, with the same caveats as above |l2j. With 
increasing D x the exponent x a ls° grows (and z de- 
creases). Obviously this cannot happen indefinitely. We 
estimate the upper limit of N in the following way: No- 
tice that the Eqs.Q and (J2J along with the prescribed 
noise correlations (i.e., equivalently the dynamic gener- 
ating functional) are of conservation law form, i.e. they 
vanish as k — > 0. Thus there is no information of any 



infrared cut off in the dynamic generating functional. 
Moreover, we know the solutions of the equations exactly 
if we drop the non-linear terms (and hence, the expo- 
nents: x — 1 — d/2, z = 2). Hence physically relevant 
quantities like the total kinetic and magnetic energies, 
J, (u(k, i)u(— k, t)) and J, (b(k, t)h(— k, t)), remain finite 
as the system size diverges, and are thus independent of 
the system size. Since the non-linear terms are of the 
conservation law form, inclusion of them cannot bring 
a system size dependence on the values of the total ki- 
netic and magnetic energies. However, if x continues 
to increase with D x at some stage these energies would 
start to depend on the system size which is unphysical 
[l4| . So we have to restrict N to values smaller than the 
maximum value for which these energy integrals are just 
system size independent: This gives N max = \{dj2 + 1). 
Note that the limits on N and N impose consistency 
conditions on the amplitudes of the measured correlation 
functions but not on the bare noise correlators. 

Antisymmetric cross-correlations stabilise the short- 
range fixed point with respect to perturbations by long- 
range noise with correlations cx k~ y , y > 0. This can eas- 
ily be seen: In presence of noise correlations sufficiently 
singular in the infra-red limit, i.e. large enough y, the dy- 
namic exponent is exactly given by 0,0] z\ T = — |. 
The short range fixed point remains stable as long as 
z sr < z\ t which gives y < — 2 + (l + N)d/2. Hence we con- 
clude that in presence of antisymmetric cross-correlations 
long range noise must be more singular for the short 
range noise fixed point to loose its stability or in other 
words, antisymmetric crosscorrelations increases the sta- 
bility of the of the short range noise fixed point with 
respect to perturbations from long range noises. 

We have seen that the amplitudes of the cross- 
correlation function play a quite crucial role in determin- 
ing the long wavelength properties of the system. In our 
analysis we used only short range noise, which is enough 
to elucidate the basic points. However, a Langevin de- 
scription of many systems often requires a noise term 
with correlations becoming singular in the long wave- 
length limit, such as fully developed MHD These 
systems are typically characterized by a set of anoma- 
lous exponents for higher order correlation functions. Be- 
low we give an illustrative example to highlight the ef- 
fects of symmetries on the anomalous scaling exponents 
of higher order correlation functions in the passive vec- 
tor limit where the velocity field u is assumed to obey a 
Gaussian distribution (instead of Eq.Q) with a variance 
(Ui(k,t)uj(-]si,0)) = (k2+ lf 2^)2+^/2 [atPij + Qij] where 
< e < 2, which makes the model analytically tractable. 
As before, the magnetic field b is governed by Eq.J^J. 
The tensor Pij is the transverse projection operator, Qij 
is the longitudinal projection operator. The parameter 
a > determines the extent of incompressibility of the 
u field. Thus in this problem a appears as a tuning 
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parameter in the multiplicative noise, very much like N 
and N appeared in our previous results. By following a 
field-theoretic dynamic renormalization group procedure 
in conjunction with Operator Product Expansion |l6| . we 
calculate the scaling exponents of the structure functions 



S n (r) = ([#(x + r)-0(x)] 2 ") ~ 
one-loop approximation we find 

Cn = 2n 



<" C ™, b = V0. Within a 



da + 1 — a 



. a 3(1 — a) 



(16) 



This clearly demonstrates that even for the linear prob- 
lem there is an continuous dependence of the scaling ex- 
ponents on the parameter a, characterizing the extent to 
which the velocity field is compressible. We expect this 
to hold also for the nonlinear problem, whose analysis is 
significantly more complicated. 

Let us now review our results in the context of 
some some physically relevant systems. Our results 
are relevant for a wide class of non-equilibrium sys- 
tems. In MHD turbulence the cross-correlation function 
(tii(k, t)bj(— k, t)) is, in general, nonzero |l7j . and as be- 
fore, is odd and imaginary in k. Similar calculations as 
here for MHD show that two dimensionless numbers, the 
magnetic Prandtl number P m and the ratio of the mag- 
netic to the kinetic energy, are non-universal; they are 
functions of both the symmetric and antisymmetric part 
of the cross-correlation amplitudes. Another system of 
interest is the dynamics of a drifting polymer through a 
solution (lif . Here the hydrodynamic degrees of freedom 
are the transverse and longitudinal displacements with 
respect to the mean position. Dynamic light scattering 
experiments can be preformed to investigate the effects of 
cross-correlations discussed here. Our results are signif- 
icant also for coupled growth of nonequilibrium surfaces 
|2cj . and sedimenting lattices 0,0]. 

In summary, we have demonstrated that cross- 
correlations between two vector fields can drastically al- 
ter their asymptotic statics and dynamics at long length 
and time scales. The symmetric and anti-symmetric part 
of the noise cross-correlation function have different ef- 
fects. The symmetric part leaves the scaling exponents 
unaffected but yields amplitude ratios of the various cor- 
relation functions, which continuously depend on the am- 
plitude of the noise cross-correlation (see Eq.(12)). In 
contrast, the asymmetric part leaves the amplitude ra- 
tios unaffected, but leads to continuously varying expo- 
nents (see Eq.(15)). In both cases the continuos vari- 
ation with the noise amplitude of the cross-correlations 
is not arbitray but constrained by scaling relations (see 
Eqs.(12) and (15)), a feature, present also in our results 
on the multiplicative noise driven linear system. We 
have shown this using renormalization group methods 



and mode coupling theory, confirmed by some prelimi- 
nary simulations [l^- Recently, Drossel and Kardar 0] 
have studied a set of coupled Langevin equations describ- 
ing the interplay between phase ordering dynamics in the 
bulk and roughening dynamics of the interface of binary 
films. They find a similar continuous variation of the dy- 
namical exponent with the coupling strength of the bulk 
and surface fields. Nonperturbative analysis or numeri- 
cal simulations may be necessary to resolve the questions 
about the rough phase more satisfactorily. In the light of 
our results it might also be interesting to examine the ef- 
fects of cross-correlations on the multiscaling properties 
of MHD in experiments and/or numerical simulations. 

One of the authors (AB) wishes to thank the Alexan- 
der von Humboldt Stiftung, Germany for financial sup- 
port. We thank S. Ramaswamy, J.K. Bhattacharjee and 
S.M. Bhattacharjee for useful discussions. 



[1] 
[2] 



[3] 
[4] 
[5] 
[6] 
[7] 

[8] 



[9] 

[10] 

[11] 
[12] 



[13] 
[14] 



[15] 

[16] 
[17] 



[18] 
[19] 
[20] 
[21] 



U.C. Tauber et al, Phys. Rev. Lett. 88, 045701 (2002). 

B. Schmittman and R.K.P. Zia, in Phase transitions and 

critical phenomena, eds. C. Domb and J.L. Lebowitz, Vol. 

17 (Academic Press, London, 1995). 

U.C. Tauber, Adv. in Solid State Physics 43, ?? (2003) 

U.C. Tauber and E. Frey, Europhys. Lett. 59, 655 (2002). 

D. Erta§ and M. Kardar, Phys. Rev. Lett. 69, 929 (1992). 

A. Levine et al., Phys. Rev. Lett. 81, 5944 (1998). 

R. Lahiri and S. Ramaswamy, Phys. Rev. Lett. 79, 1150 

(1997). 

S. Yanase, Phys. Plasma 4, 1010 (1997); J. Fleicher and 
PH. Diamond, Phys. Rev. E 58, R2709 (1998); A. Basu, 
J.K. Bhattacharjee and S. Ramaswamy, Eur. Phys. J. B 
9, 725 (1999). 

D. Forster, D. R. Nelson, and M. J. Stephen, Phys. Rev 
A, 16, 732 (1977). 

A. Basu, Phys. Rev. E 61, 1407 (2000). 

J.K. Bhattacharjee, J. Phys. A, 31, L93 (1998). 

E. Frey, U.C. Tauber, and H.K. Janssen, Europhys. Lett. 
47, 14 (1999); H.K. Janssen, U.C. Tauber and E. Frey, 
Eur. Phys. J B 9, 491 (1999). 

A. Basu and E. Frey, unpublished. 

In contrast with our model system, physical systems, 
such as fluid and MHD turbulence, where total energy 
diverges with the system size are described by stochas- 
tically driven Langevin equations with noise correlators 
diverging with the system size. 

A.K. Chattopadhyay and J.K. Bhattacharjee, Europhys. 
Lett., 42, 7381 (1998). 

L.T. Adzhemyan et al, Phys. Rev. E, 58, 1823 (1998). 

A. Basu, A. Sain, S.K. Dhar and R. Pandit, Phys. Rev. 
Lett. 19, 19 (1999); A. Basu, PhD Thesis, Indian Insti- 
tute of Science (1999). 

D. Ertas. and M. Kardar, Phys. Rev. E 48, 122 8 (1993). 

B. Drossel and M. Kardar, cond-mat/0307579 (eprint). 
AL. Barabasi, Phys. Rev. A 46, R2977 (1992). 

D. Das et al, Phys. Rev. E 64, 021402 (2001). 



